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Abstract 

The effects of higher order corrections in ultra-relativistic nuclear collisions 
are considered. It is found that higher order contributions are small at low 
energy, large at intermediate energy and small again at very high energy. An 
explanation for this effect is given. This means that the Weizsacker- Williams 
formula is a good approximation to use in calculating cross sections and beam 
lifetimes at energies relevant to RHIC and LHC 

PACS numbers: 25.70. Jj 



Heavy ion and relativistic heavy ion accelerators have so far been con- 
structed with an incident beam and a fixed target. The first heavy ion 
collider, named the Relativistic Heavy Ion Collider (RHIC) is currently 
under construction at Brookhaven National Laboratory. This machine will 
collide heavy ions all the way from the proton up to Pb nuclei with an energy 
per beam of 100 GeV/nucleon. This is equivalent to a fixed target acceler- 
ator with a beam energy of 21.7 TeV/nucleon. The Large Hadron Collider 
(LHC) is currently being planned for the purposes of proton-proton colli- 
sions. However LHC will also be capable of accelerating heavy ions. One of 
the main purposes in constructing these two relativistic heavy ion colliders 
is to undertake a search for the quark gluon plasma phase. 

Heavy nuclear beams installed in heavy ion colliders will only have a finite 
lifetime due mainly to beam-beam Coulomb dissociation and electron pair 
production. In designing heavy ion colliders it is important to have good 
estimates for the various beam lifetimes 0. Such an estimate for U-U beams 
installed in RHIC have been carried out || . However in that work the effects 
of higher order electromagnetic processes || was left out. Such effects are 
known to be small (about 4%) at low energy |3J]. However these effects have 
been found to be larger (about 12%) at higher energies ||. Thus the question 
immediately arises as to how big these effects might become at the very high 
energies planned for RHIC and LHC. In particular, if the effects are quite 
large then previous estimates of the beam lifetimes |], might be incorrect. 
This may have a significant effect on accelerator operations. 

Generally the higher order cross sections are smaller than the Weizsacker- 
Williams cross section. In multiple electromagnetic processes the excitation 
of the single photon giant dipole resonance is depleted by multiple higher 
order electromagnetic processes, which can excite double and multiple giant 
resonances, generally lying at double and multiple integral values of the en- 
ergy of the single dipole state. The Weizsacker- Williams calculation assumes 
that all of the photons excite only the single giant dipole resonance state. 

Calculations were performed for the Coulomb dissociation cross sections, 
o~wwi using the Weizsacker- Williams method |J for Uranium- Uranium (UU) 
collisions at a variety of energies. The calculations were repeated for the 
cross sections including higher order corrections, a at, as discussed in M. The 
percentage difference, 100(aww — o~n)/o-ww, is then calculated. The results 
for fission of one of the U nuclei are presented in Figure 1. As discussed 
above, it is seen that the higher order effects are small at low energy and 
then grow as the energy is increased. However the surprising feature is that 
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as the energy is increased even further, the higher order effect starts getting 
smaller again. For RHIC and LHC energies the percentage difference is less 
than 2%. This is good news for calculations of beam lifetimes, because it 
means that the old Weizsacker- Williams calculations |2| remain correct. 
It has been found previously || that other corrections, such as Rutherford 
bending and electric quadrupole effects, can also be neglected at high energy. 
Thus at the very high energies of RHIC and LHC one can essentially ignore 
all corrections (including higher order corrections), and just rely on the basic 
Weizsacker- Williams formula. The calculations represented in Figure 1 are 
for fission, but similar results were also found for other Coulomb dissociation 
reactions such as single and double nucleon removal and also for other nuclei. 

Further analysis revealed the source of the behavior seen in Figure 1. In 
performing a Weizsacker- Williams calculation one folds and integrates the 
virtual photon spectrum with the photonuclear cross section as in 

<?ww = J N(E)a(E)dE (1) 

where a is the nucleus- nucleus Coulomb dissociation cross section, N(E) is 
the virtual photon spectrum as a function of virtual photon energy E, and 
<j(E) is the photonuclear cross section. The virtual photon spectrum is a 
smoothly dropping function of virtual photon energy || , while the photonu- 
clear cross section is a bell shaped curve centered on a fixed photon energy. 
As the heavy ion energy is increased the amount of "overlap" between the vir- 
tual photon spectrum and the photonuclear cross section steadily increases. 
In impact parameter space this cross section is written |J 

a = 2tt J P(b)bdb (2) 

where P(b) = PWW(b) is the Weizsacker- Williams probability given by 

PWW{b) = $(&) (3) 

where 

$(b) = J n(E,b)a(E)dE/E (4) 

with n(E, b) being the impact parameter dependent virtual photon spectrum 
||. Thus if PWW(b) is used for P{b) in equation (fj) then one reproduces 
the numerical values obtained from (|l|). The higher order cross section is 
calculated with equation (|2]) except that 

PN(b) = $(6)e-* (6) (5) 
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is used for the probability P(b). This Poisson probability distribution results 
from the statistical independence of the emission of a number N of photons 
emitted from a classical current source 

Now to explain Figure 1. The probabilities, multiplied by b, (since it is 
b x probability that gets integrated in (H)) are plotted in Figures 2-4. At 
very low energies <J>(6) will be very small, so that $(&) and $(6)e -<1>( '^ will be 
approximately the same leading to equality of the Weizsacker- Williams cross 
section o~ww and the higher order cross section <tat. This is seen in Figure 
2, and is reflected in Figure 1 at the low energy end where the Weizsacker- 
Williams and higher order cross sections are approximately the same. As 
the energy is increased however, $(6) increases and $(&)e~*( b ) starts getting 
smaller than $(6), as shown in Figure 3, meaning that is smaller than 
o~ww which explains the rising part of Figure 1 between 1CT 2 GeV/nucleon 
and 10° GeV/nucleon. Remember that it is the areas under the curves in 
Figures 2-4 that gives the cross sections according to (H). 

Figure 4 shows the probabilities (times b) at very high projectile energy. 
Even though there is still a difference between the probability curves, however 
now both curves extend out to much larger impact paramteres where the 
differences are much smaller. (High energies probe larger impact parameters.) 
This results in the total areas under the two curves in Figure 4 being much 
closer in value to each other than the two curves in Figure 3. (This was 
seen by explicit computation.) This means that for the case of Figure 4 the 
resulting cross sections (calculated from (Q)) are in fact quite close to each 
other. This explains the small percentage difference value in Figure 1 at high 
energy. Consequently, the curve in Figure 1 has reached some maximum 
value (near 1 GeV/nucleon) and then starts falling as the energy increases. 
These various effects then explain the behaviour seen in Figure 1. 

At very high energy then, one can safely use just the Weizsacker- Williams 
formula for calculating beam lifetimes, even for the heaviest of nuclei where 
one might have thought that higher order effects might have been important. 
As an application of this a Uranium beam lifetime is calculated for LHC 
The method of calculation is the same as described previously 0. The 
beam parameters || are listed in Table 1, where k is the number of beam 
intersections, Lq is the initial luminosity, B is the number of bunches and 
Nb is the number of particles per bunch. Cross sections and lifetime results 
are listed in Table 2, where some previous results || || for RHIC are also 
listed for comparison. (The greatest uncertainty in these calculations is the 
pair production cross section a e + e - ||. Also in @, the value for a e + e - should 
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have been 85 barn.) Table 2 shows that Uranium beams will live sufficiently 
long so that they could be installed at RHIC and LHC. 

In conclusion, it has been shown that higher order contributions to Coulomb 
dissociation are small at low energy, large at intermediate energy and then 
small again at very high energy Consequently, the use of the Weizsacker- 
Williams formula in calculating cross sections and beam lifetimes is expected 
to be accurate at RHIC and LHC energies. 

This work was supported by NASA grant 97-367-NM. 
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Beam parameter 


RHIC 


LHC 


k 

L (cm' 2 s~ r ) 
B 

N B 

energy/nucleon/beam (GeV) 


6 

9.2 x 10 26 

57 
1.1 x 10 9 

100 


1 

2 x 10 27 

608 
9.4 x 10 7 
2760 



Table 1: Beam Parameters 







a (barns) 






t (hours) 




fission 


IN 


2N 


e + e 




RHIC 


33 


48 


31 


85 


11 


LHC 


56 


80 


52 


135 


17 



Table 2: Weizsacker- Williams Results 



6 



1 0" 3 10" 2 10" 1 10° 10 1 10 2 10 3 10 4 

Tlab (GeV/nucleon) 



Figure 1 

Percentage difference between Weizsacker- Williams and higher order cross 
sections for fission in UU collisions as a function of projectile kinetic en- 
ergy for an equivalent fixed target accelerator. The RHIC energy (21,700 
GeV/nucleon) and the LHC energy (16,237 TeV/nucleon) lie beyond the far 
right of the graph. 
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Figure 2 

Impact parameter multiplied by probability as a function of impact parameter 
for Weizsacker- Williams and higher order theory for a projectile energy of 
50 MeV/nucleon. The area under each curve gives the total cross section. 



8 



3 




20 25 30 35 40 45 50 
Impact Parameter, b (fm) 



Figure 3 

Same as Figure 2 except for a projectile energy of 1 GeV/nucleon. 
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Figure 4 

Same as Figure 2 except for a projectile energy of 1 TeV/nucleon. 
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